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ABSTRACT 

The effect of solar wind on dust particle, meteoroid, is investigated. Rotation of the 
particle is also considered. Detail derivations of both equation of motion and Euler's 
dynamical equations, are presented. The most simple form of the action of the solar 
wind is considered and the derived equation of motion reduces to "pseudo-Poynting- 
Robertson effect" for non-rotating spherical particle. The final equation of motion 
contains also angular velocity term and it may cause change of orbital plane of the 
meteoroid. Analysis of the Euler's dynamical equations shows that rotational axis of 
the meteoroid changes in time and equilibrium state for the orientation of the axis 
of rotation does not exist for spherical particle. 
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1. Introduction 

At the beginning of the past century, Poynting (1903) has formulated a problem of finding 
equation of motion for a perfectly absorbing spherical particle under the action of electro- 
magnetic radiation. Robertson (1937) derived a correct equation of motion for a perfectly 
absorbing spherical particle. This result has been applied to astronomical problems for 
several decades and is known as the Poynting- Robertson (P-R) effect. Relativistically 
covariant equation of motion for general properties of spherical particles was presented 
by Klacka (2000, 2004). Complete understanding of the action of the incoming electro- 
magnetic radiation on the motion of the spherical dust particle was obtained by Klacka 
(2008a, 2008b) and Klacka et al. (2009a). Partly similar effect on interplanetary dust 
particles is caused not only by solar electromagnetic radiation, but also by stream of 
charged particles flowing from the Sun. This stream continuously emanating from the 
Sun's surface consists mostly from electrons and protons, a— particles and heavier atoms. 
The stream is called also corpuscular solar radiation, or, shortly, solar wind. Dynamical 
effects of solar wind on motion of meteoroids are similar to dynamical effects of solar elec- 
tromagnetic radiation and they are known also as "pseudo-Poynting-Robertson effect" 
(Dohnanyi 1978). Detail understanding of the effect of solar wind on motion of spherical 
interplanetary dust particle is presented in Klacka et al. (2009b). 

All the mentioned equations of motion for dust particles under action of electromag- 
netic radiation and solar wind are correct for non-rotating dust grains. What is the effect 
of rotation? How does dust particle rotate? In order to be able to answer the questions, 
we will deal with the effect of rotation on the dust particle. We will obtain results not 
only for time evolution of rotation of the particle, but also equation of motion of the 
rotating particle. 



2. Proper reference frame of the dust particle - stationary particle 

We will use considerations and accuracy to the first order in v/u and ljR/u (v is helio- 
centric velocity of the spherical dust particle of radius R, u is the speed of solar wind 
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with respect to the Sun and us is angular velocity of the rotation of the particle about a 
rotational axis containing center of mass of the particle). 

The term "stationary particle" will denote a particle which does not exhibit transla- 
tional motion in a given inertial frame of reference. Primed quantities will denote quan- 
tities measured in the proper reference frame of the center of mass of the particle. The 
particle may be characterized with a rotational motion in the proper reference frame of 
the center of mass. 

Let S' is an incident flux density of solar wind energy (energy flow through unit area 
perpendicular to the ray per unit time) and geometrical cross section of the spherical 
particle is 7r R 2 . We have S' = n' mi \u'\ 2 c, where n' is concentration (number of particles 
per unit volume) of solar wind particles of mass mi (for simplicity of these considerations 
only one type of particles is considered) , c is the speed of light, u' — u' S' is velocity vector 
of solar wind particles: unit vector S' denotes direction and orientation of the incident 
beam of solar wind particles measured in the proper reference frame of the center of 
mass of the dust particle. Relations (1) - (19) in Klacka and Saniga (1993) are significant 
for the effect of incident solar wind particle. Moreover, Lorentz transformation yields for 
elementary force acting on the elementary area of the spherical particle due to action 
of the incident solar wind (dp' / dt) el . inc = {dp" / dt) el . inc + E'l nc v' rot .J c 2 , where double 
primed quantities are measured in the frame of reference connected with the rotating 
(rotational velocity v' rot . el ) elementary area on the sphere. We have (dp" / dt) el . inc — n' el 
[(dA') e i cos a'] \u' el \ mi u' el , n' el = n' = n (within the considered accuracy), u' el = u' — 
w x If we take all these results into account, then we can write for an elementary 
force caused by the incident solar wind 

] = n' mi u' 2 [{dA') el cos a'] x 



d t 



el;inc 



\U'\ J \w 

+ ri mi \u'\ [{dA') el cos a'] oj x £' , (1) 

where is position vector of the clement on the surface of the spherical particle with 
respect to the center of mass of the particle, (dA') el is area of the element. Within the 
considered accuracy, Eq. (1) is equivalent to 



fdp^\ = S' (dA') el cosa' ( gl \ g , 

\dtj e , inc c { \u'\ J 



(2) 



if the incident flux density of solar wind energy S' — n' mi \u'\ 2 c is used: this enable 
to consider more general composition of solar wind, as for various types of solar wind 
particles. 

Eq. (2) enables to find the total force of the incident solar wind on the dust particle, 
and, also the total torque acting on the particle: 

" b / total;inc e ; \ u L / el;inc 

M' tota , inc = -£t'x(^j) . (4) 



In order to find (dp' /dt) total . inc and M' total . inc , we may perform calculations in some 
special frame of reference: 

(dA') el = [(R sin a') d X'] R d a' , 
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£' = R (cos a', cos A' sin a', sin A' sin a') , 
S' = - (1, 0, 0) , 

OJ = (U> x , UJy, u> z ) , 

a' 6 (0, 7T / 2) , 

A' e (0, 2 tt) . (5) 
Putting Eqs. (2) and (5) into Eqs. (3) and (4), using also ^2 el — ► J, one finally obtains 

^) = ^ 5' , (6) 

d * / totaHnc C 

M 'totamnc = -\ {« " (« ' «'} • (7) 

As for the "outgoing part" of solar wind, we will suppose that the solar wind does 
not impart any momentum to the corresponding element of the area on the surface of 
the dust particle. Thus, we have (dp" /dt) el . out — for elementary force acting on the 
elementary area of the spherical particle due to action of the "outgoing" solar wind. 
Using the relation (dp' / dt) el . out = (dp" / dt) el . out + K ut v' rot .J c 2 , we can write 

dp'\ x' S' u> x £' 

[(dA') el cos a'} — -f- , (8) 



dt Li-.out c ei |u 



where it is supposed that the x' — th part of the incident energy per unit time is lost 
from the dust particle (see Eq. 20 in Klacka and Saniga 1993). 

Eq. (8) enables to find the total force acting on the particle due to the "outgoing" 
solar wind, and, also the total torque acting on the particle (the sign minus is due to the 
conservation of momentum): 




dp' 
~~dt 



(9) 



ehout 



M' total;out = ^'x\ 1 • (10) 



, , , ^ t / ehout j 

Detail calculations in a special frame of reference, e. g. given by relations (5), yield 
dp'\ x' S' tt R 2 2 R 

u> x S , (11) 



dt J total-out C 3 I"'. 

M' tota ,, out = -\ X ' S c '** {3 u, - ( W • S') S'} . (12) 
The total force and torque acting on the particle due to solar wind: 

d ^ > d ^ + m\ , (i 3) 



d t J total \ dt J tota i. inc \ dt / total . out 



M' tota i — M' total . inc + M' total . out , (14) 
or, using Eqs. (6)-(7) and (11)-(12): 

dp'\ S' 7r R 2 ( f , 2 R _,1 

1 S' + x' - — u x S' } , (15) 



d t 



total 



M 'totai = ~\ {(1 + 3 x') « - (1 + x') (u, • S') S'} . (16) 
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3. Stationary frame of reference - equation of motion 

By the term "stationary frame of reference" (laboratory frame) we mean a frame of 
reference in which particle moves with a velocity vector v = v(t). The physical quantities 
measured in the stationary frame of reference will be denoted by unprimed symbols. 

Our aim is to derive equation of motion for the particle in the stationary frame of 
reference. We will use the fact that we know this equation in the proper frame of reference: 
Eq. (15). In order to find the equation in the stationary frame of reference, we will use 
the Lorentz transformation for force acting on the dust grain: 

dp = dp' E' inc - E' out v 

dt dt c c ' [ ! 

Since 

E' - t' E' 

'-'out •"' '-'vac 1 

E' mc = n' mi | it' | c 2 tt R 2 , (18) 

(compare with the right-hand side of Eq. 1: E' inc = J2 e i n ' m i \ u '\ ^ \{dA') e i cos a']), we 
have 

< LP = ^+ (1 - x') ri mi \u'\ tt R 2 v . (19) 

Since 5" = S (1 - 2v ■ S/\u\), S' = (1 + v ■ S/\u\) S - v/\u\, Eqs. (15) and (19) yield 

dp Stt R 2 f / v -5\ v , 2 R a \ 

dt c { \ \u\ J \u\ 3 \u\ J 

,. S tt R 2 v . . 

+ (1 - x } R ' (20) 

The last term in Eq. (20) corresponds to the change of grain's mass, in accordance with 
Eq. (18): 

^ = (i - S -^r ■ (2i) 

dt v ' c\u\ K ' 

Since dp/dt = d(mv)/dt — m dv/dt + v dm/dt, Eqs. (20)-(21) lead to 
dv S tt R 2 f/ v-S\ „ v ,2R 

di mc \V \u\ J \u\ 3 \u\ J v ; 

Mass of the homogeneous spherical particle is m = ( 47r/3 ) g i? 3 , where g is mass 
density. Eq. (21) reduces to 

^ - - H (23) 
where r is distance between the particle and the Sun, and 

K={x' - 1) / 2 ^ . (24) 
4 g c \u\ 

is a constant for a given material. According to Dohnanyi (1978) and Kapisinsky (1984) 
we can take K w 4 x 10~ 9 cm year" 1 (see also Leinert and Grun 1991). The phenomenon 
presented in Eq. (23) is known as the "corpuscular sputtering". The units used in Eq. 
(23) are: [ R ] = cm, [ t ] = year, [ r ] = AU, [ K ] = cm year^ 1 . Eq. (24) enables to 
find the vale of the dimensionless factor x'\ x' — 1.9 x 10 11 K[cm year' 1 ] g[g cm' 3 ] + 
1. As an example, for the case K = 4 x 1CT 9 cm year' 1 and g = 1 g cm' 3 , we obtain 
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4. Rotational motion of spherical particle 

On the basis of Eq. (16), we can immediately write Euler's dynamical equations: 
d [(2/5) m R 2 tJ\ 1 Sir R 4 



dt 4 c \u\ 



{(1 + 3i')«- (1 + x 1 ) (w • S') S'} , (25) 



where we have used S instead of S' and \u\ istead of \u'\, within the considered accuracy. 
The moment of inertia of a sphere of uniform density is (2/5) m R 2 . 

In order to solve Eq. (25), together with Eqs. (22)- (23), we use, as standardly, Euler's 
angles -0, "& and <p>. We have 

d ojy 5 ir S R 2 „ ,. . ... _ , 5 dR 

^ = -^^{(i + 3*') ^ - (iW) (w .fir)M - |^ 2 <, 

at 8 mc|u| R dt 

d u>3' 5 7r S R 2 „ ,. , ,. . _ , 5 dR ,„„. 

-7T = - X {(1 + 3x) W3 ' - + ( "' 5) * 3 ' } - * 77 ^ ' (26) 

where we have used the fact that scalar product satisfies Ylp=i w p<Sp — Sq=i w q >S q ', and, 
S 1 = Sii + + in the stationary frame of reference, S' = Syi' + Syj' + Syk' 
in the frame of reference rotating with the particle. Moreover, the well-known relations 
are: 

u)v = tp simp simd + z9 cos<p , 
oj2> = tp cos <p sin 1? — $ sm<p , 

U31 — tp COS I? + (p , (27) 

Ui = <p simp sin 1? + $ cosip , 
ui 2 = — <p cos?/; sint? + $ simp , 

U3 = tp cost} + tp , (28) 

S = Si i + S 2 j + S3 k , 

S = -* + -j + -fe, 

r r r 

r = 11 + y j + zfc, (29) 

where r is position vector of the particle with respect to the Sun, source of solar wind, 
r = \r\, v = dr/dt, and, 

(Sv \ ( Ru R12 R13 \ ( S\ \ 
S*> = i?21 #22 i?23 ^2 , (30) 

-5*3' / V #31 #32 #33 / \ S3 / 

Rn — cosip cos ip — cosfl simp simp , 

i?i2 = sintp cos ip + cos 1? cosip simp , 

R13 = sin-# sin ip , 

#21 = — cosip sintp — cos $ simp cos ip , 

R22 = — simp simp + cosfl cos?/; cos<p , 

R23 — sin 1? cosip , 

R31 = sin 1? simp , 

#32 = — sim) cosip , 

#33 = cost? . (31) 
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Eqs. (28) and (29) yield 

u> ■ S = (j} cos ip + if sin i9 sin tp^j — + (^9 sin ip — (p sin z9 cos ipj — + 

+ (ip cos?? + - . (32) 
We can take into account a non-dimensional parameter, for Solar System 



r 2 Su gh t tt R 2 -, _ L Q tt R 2 -, _ OA ^, n - 3 R 2 [ 
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= ' _ J " q' = ^ " " Q' = 2.4 x 10- 3 Q' , (33) 

where Lq is the rate of energy outflow from the Sun, the solar luminosity, Mq is 
mass of the Sun and G is the gravitational constant and Q' pr is a dimensionless ef- 
ficiency factor for radiation pressure (Klacka 2004). The non-dimensional parameter 
("the ratio of radiation pressure force to the gravitational force") reduces to (3 = 
5.7 x 10~ 5 Qpr / (g[g / cm 3 ] R[cm]), for homogeneous spherical particle: g is mass den- 
sity and R is radius of the sphere. On the basis of Eq. (33), we can immediately write 
for the quantities of solar wind: S tt R 2 j (m c) = f3 (r) / Q' pr ) (G Mq / r 2 ) \u\/c, where 
J7 ~ 1/3. Inserting this result into Eq. (26), using also Eq. (24), one immediately obtains 



V 


GM Q I 




+ 11 


Qpr 


r 2 ) 




3 


_n_ 


GM Q 1 


u> 


+ 11 


Q'pr 


r 2 ) 




3 


V 


GM Q j 


fx' 


+ 11 



(l + x') (u-S) Si, | , 



d UJy 5 

d t ~ 8^? 

d o> 2 ' 5 
d t ~ 8^ 

- -I't^ i^" 3 ' - {1 + x ' } ( "' S) Sa 'J ■ (34) 

The orbital and rotational motion of spherical dust particle under action of solar wind 
and gravity of the Sun is given by Eqs. (23), (27)- (34), together with equation 

£__£*a s+/ , » £Lj<2{(M _ ^) S -? + ,^„ xS l (35 ) 

dt r z Q' pr r z [ \ c c J c 3c J 

Equation of orbital motion represented by Eq. (35) depends on angular rotational 
velocity of the particle. 



5. Conclusion 

The action of solar wind on motion of arbitrarily shaped dust particle would produce 
planar motion for non-rotating meteoroid orbiting the Sun. If angular velocity of rotation 
of the meteoroid is large enough, then rotational term may change orbital plane (see Eq. 
35 for spherical particle; see also Eq. 8). Euler's dynamical equations show that orbital 
axis of the particle changes and no equilibrium for the orientation of the axis exists. 

Equation of motion and quantities for Euler's dynamical equations for arbitrarily 
shaped dust grain are given by Eqs. (2)-(4), (8)-(10), (13)-(14), (17) and the first part of 
Eq. (18). 
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